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Abstract: In this study, the use of taylor series method in the calculation of some means with single auxiliary variable 
developed in a simple random sampling and mean square error unit ratio estimators having certain properties was 
investigated. An application was performed in respect thereof. The study population (mass) included 111 secondary schools 
from 18 districts of Trabzon province. Auxiliary variable (x) was taken as the number of students whereas the main variable 
(y) was taken as the number of teachers. Sample size was calculated as 45 for unit ratio that has certain features. Afterwards, 
theoretically proposed mean and units ratio estimators having certain properties were compared numerically. Random 
sampling was performed using the SPSS 20 program thus giving an equal chance to the units sampled and variability in the 
population was protected. 

Keywords: Simple Random Sampling, Taylor Series Approach, Mean, Proportional Estimators, Mean Square Error, 

Units Ratio having Certain Features, Ancillary Information 



1. Introduction 

Proportional estimator method is not a sampling method 
but a method for finding estimator. In the case where the 
correlation between two variables is positive, information 
from auxilary variables is used. Newly proposed ratio 
estimators are frequently encountered in the literature. Data 
related to auxiliary variable were examined in detail by 
Cochran (1977). Singh-Tailor (2003) estimation is based on 
mean ratio estimation proposed in the case where mass 
correlation coefficient (p) of main and auxiliary variables 
is known. Kadilar and £ingi (2004) proposed new ratio 
estimators by implementing linear regression estimator 
rather than simple estimates and using coefficient of 
variation information (C*) from auxilary variable. Kadilar 
and Cingi (2006) developed a new mean ratio estimator by 
typing regression estimation itself in Singh and Taylor 
estimator and using coefficient of mass variable from 
coefficient of mass correlation and auxiliary variable rather 
than simple estimator. Naik-Gupta estimator (1996) is the 
unit ratio estimator having a certain property. Singh et al. 
(2008) proposed an estimation method based on the 



implementation of linear regression estimation having 
certain features rather than simple conventional estimation 
in Naik-Gupta estimation and the use of data such as 
variation coefficient and correlation coefficient. 

In our application, the estimators proposed in the 
literature were compared numerically. 

2. Taylor Series Approach 

Observation of non-linear estimates is frequently 
encountered in sampling. It may not be possible to calculate 
the variance of these non-linear estimators exactly. 
Estimations such as ratio estimate or estimation of 
regression coefficient may be used. To calculate the 
variance of a non-linear estimator, observations of non- 
linear estimator are approximated to a linear function and 
an appropriate variance equation is written. This 
linearization is usually performed by Taylor Series Method. 
In other words, Taylor series method can be used in the 
calculation of variance of an estimator. 
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In general, it can be given for the variable k as follows: 

hiY^.-.Yk) = A(li. V ,Y k ) +Y.J = 1 dj( 9 ] - Yj) + R k (Y k , a) ( 1 ) 

where h(Y 1 , Y 2 ■■■ , Y k ) is a non-linear function and h(Y 1 , Y 2 : 
■■,Yk is the estimation of this function. 



dj — 



dh(a^a 2 ,-,ci k ) 



da; 



( 2 ) 



Yl ’ Y 2h(? 1 ,?2.~,? k );" Yk 

where d represents the notation of partial derivative and 

R*(Y k ,a) = ~ Yj)(? t - Y t ) + O k 

(3) 



y 0 = I X = RX (9) 

Here it is equal to R = - and y is the sample mean of the 
variable of interest and x is the sample mean of the 
auxiliary variable. We assume that mass mean of auxiliary 
variable x is known, 

x = li = i x t ve y = Zp =1 y t (10) 

where n is the number of sample units. 

Mean square error of this proposed estimator was 
obtained as 

MSE{y 0 ) a A(S y 2 - 2 Rp yx S y S x + R 2 S X 2 ) = X (s y 2 + R 2 S X 2 ( 1 - 20 )) ( 11 ) 



Here, (^represents the terms after second order in Taylor 
series method. 

Here, the term R k (Y k , a) does not change the value of 
mean square error dramatically; but R k (Y k , a) is important 
to calculate the bias. In simple random sampling, second 
order terms, in other words, the term R k (Y k , a) is 
frequently neglected when the sample size is large enough. 
Therefore, the mean square error can be obtained by Taylor 
Series Method as follows; 



r 

Where = p — , C y and C x are the variable of interest and 

Cx y 

coefficient of variation of the auxiliary variable respectively. 

l-f 

A = — is the sampling fraction f and it is obtained by 
division of sample width to the mass width. In other words, 

n Y 

= - ; is the number of N mass units. R = - is the mass 

N X 

ratio; S y 2 and S x 2 are the mass variance of the variable of 
interest and auxiliary variable respectively [2], 



h{Y Y 2 , , Y k ) = h(Y lt Y 2 , ■■■, Y k ) + Zf =1 dj(?j - Y } ) (4) 

By using Equation (4), mean square error is obtained as 
follows; 

E[h&, Y 2 ,-,Y k ) - h(Y ll Y 2 ,..Y k )] 2 S£ E[l!f =1 dj(9j - Yj)] 2 (5) 
MSElhCY^.-.Yk)] = E[j* =1 dj(?j - Yj)] 2 (6) 



3.3. Investigation of Singh and Tailor Ratio Estimation 
and Mean Square Error 

If the mass correlation coefficient of variables of interest 
and auxiliary variable (p) is known, ratio estimate 
proposed by Singh-Tailor can be defined as follows; 

y S T = ^-(x + p) (i2) 



3. Mean Square Error of the Proposed 
Estimators 



Mean square error of this estimator is expressed by 

MSE(y ST ) = A (5 y 2 + R 2 S x 2 w(w - 20)) (13) 



3.1. Simple Conventional Estimation and Mean Square 
Error 



where w = ve 6 — — 

X+p c x 



[3], 



When n-width sample is taken from n-width mass via 
simple random sampling method, variance of y can be 
obtained as follows; Variance of mean estimator, 

ny) = £(y-F) 2 =M^l (7) 

Where f = ^ ve A = ^ is taken, Equation (8) becomes 

ny) = ^^ = ^ 2 y = As y 2 (8) 

Sy 2 = Zg=1 ^ y)2 [11 

3.2. Conventional Propotional Estimation and Mean 
Square Error 

Conventional ratio estimation for ( Y) mean mass of 
variable of interest y is 



3.4. Some Estimators Proposed by Kadilar-Qngi and 
Mean Square Error 

Kadilar and (^mg 1 (2004) proposed the following ratio 
estimators. 

5W = (14) 

y K p = ^f^[x + c x ] (15) 

R K pi = y+b(x x) and b = s x 2 represents the sample 
variance of auxiliary variable and s yx is the sample 
covariance between variables of interest y and auxiliary 
variable x. General form of this expression was proposed 
by Ray and Singh (1981) [4]. Proportional estimator 
proposed by Ray and Singh is as follows; 



- _ y+6(x“-x“) ^ 

Y r s ~ “T A 



( 16 ) 
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In this estimator, if b = 0 and ^ = 1 , it becomes 
conventional proportional estimation. In addition, in the 
case of a = 1 and ^ = 1, proportional estimator proposed 
by Ray and Singh corresponds to that proposed by Kadtlar 
and Qngi. Mean square error of Kadtlar and (,'ingi 
estimator is expressed as follows; 

MSE{y K Ci ) = X[(B + V + S y 2 - 2 (B + R K(;i )S xy )] (1 

7) 

Here, if B = and is squared 

It is expressed as 

MSE(y KV ) = A[(R kq 2 S x 2 + 5 y 2 (l - p 2 )] ;l = 1,3 (18) 

where R KQ1 = R = j , R K(3 = [5] 

Inference: In Equation (17), if B = — R K p, mean square 
error of y K ^i estimator is equal to the y variance of simple 
estimator. 

If B = — R K?i , 

MSE(y KCi ) = AS y 2 = V(y) 

Kadtlar and £tngt (2006) developed a new ratio 
estimator by utilizing regression estimation rather than 
simple estimation in Singh and Taylor Estimator and using 
coefficient of mass correlation and variation coefficient of 
auxiliary variable. This estimator is as follows; 



y+b(X-x) y 

yprl= *- +p (* + P) 


(19) 


y+b(X-x) . y 

y Vr2~ xCx + p ( XC x + P) 


(20) 


y + b(X-x ) „ y 

y pr 3 - xp + C x ( Xp + C J 


(21) 


Mean square error of this estimator can 


be expressed as 



follows 



A 3 

By using this Equation, P = — and p = - mass and 
sample ratio estimations and ratio equation of the desired 
events can be obtained. Here, = <^>1 and a = XP=i4 ) i 

show the total number of mass and sampling units 
respectively [1]. 

3.5.1. Investigation ofNaik and Gupta ratio Estimators 
and Mean Square Error 

Vng — p P (24) 

where y is the sample mean of the variable of interest, p 
is the unit ratio having certain features in the sampling. 
Mean square error of this estimation is 

MSE(y NG ) = A(R NG 2 S ( j ) 2 — 2R NG S y( p + S y 2 ) 
where S 0 2 = ^Zi=i(.4>i - P Y and 

S y0 =^ I I]f=i(0i-7 5 )(yt-F) [7]. 

3.5.2. The Proposed Singh, Chaucer, Sawan and 
Smarandache Estimates and Mean Square Error 
Error 

Estimations following Ray and Singh Ratio Estimations 
are 



t _ y+b^p-p) p 
1 V 


(25) 


t2 -^ [P + cA 


(26) 




(27) 


In Equation (25), if b^ = 0 , the proposed estimator 


becomes Naik and Gupta estimation in 
Mean square error of this estimator is 


Equation (24). 


MSE{t{) = AlRfS^ 2 +S y 2 (l -p 0 2 )] 


; i = 1, 2, 3 



MSE{y vri ) = A[R pr 2 S x 2 + S y 2 (l - P 2 )] /* = 1,2,3 (2 



where R prl - +p ,R pr2 ^ +p , R pr3 KpyCx [6]. 

3.5. Some Ratio Estimators that have Certain Feature 
which were Developed in Simple Random Sampling 
Using Auxiliary Information and Mean Square Error 

Here, R in ratio estimation is replaced by B in regression 
estimation. The use of auxiliary variables can improve the 
efficiency of estimation in the case of where auxiliary 
variables and the variable of interest are highly correlated. 
Suppose that n-width sample is taken from n-width mass 
(population) via simple random sampling without 
replacement. y- t and 4>i represent the observation values in 
the variables for i = 1, 2, , N respectively. Suppose that 

_ fl, if the desire feauture is achieved 
® >l l 0, other cases ' 



where S 0 2 = — Zi=i(4>i - P Y 

Sy^rhlZMi-PXyi-Y) [ 8 ], 



4. Theoretical Comparisons 

4.1. Comparison of Singh and Tailor Estimator and 
Conventional Proportional Estimation 

Conventional proportional estimator; 

MSE{y 0 ) a A(R 2 S X 2 - 2Rp yx S y S x + S y 2 ) = X (s y 2 + ff 2 S* 2 (l - 20)) (28) 

We have seen conventional proportional estimator in 
Equation (28). Conventional proportional estimation is 
obtained by inequality of mean square error of Singh and 
Tailor Estimator 




144 



Tolga Zaman et al. : Investigation of Some Estimators Via Taylor Series Approach and an Application 



A(S y 2 + R 2 S x 2 w(w - 2 6 »)) < A(S y 2 + R 2 S X 2 ( 1 - 20 )) 
w(w — 20) <1 — 20 
w 2 - 1 - 2w0 + 20 < 0 
(w — l)(w + 1) — 20(w — 1) < 0 
(w - l)(w + 1 - 20) < 0 
This inequality is satisfied when 

w — 1 < 0 ve w + 1 — 20 > 0 (29) 

or 

w — 1 > 0 vg wTl — 20<O (30) 

Under these conditions, mean square error of Singh and 
Taylor Estimator is less than that of conventional 
proportional estimator. Therefore, it is more effective. 

4.2. Comparison of Singh and Tailor Estimator and 
Simple Conventional Estimation 

A (s y 2 + R 2 S x 2 w(w - 20 )) < AS y 2 

R 2 S x 2 w 2 < R 2 S x 2 w29 

J<2 (31) 

Under condition (31) Singh-Tailor estimation is more 
effective than simple conventional estimation. 

4.3. Comparison of Kadilar-Qngi (2004) Estimation and 
Conventional Proportional Estimation 



F(y) = ASy 2 


(34) 


V(y)-MSE(ti)>0 ;i= 1,2,3 


(35) 


ASy 2 > A[R 2 S^ 2 + S y 2 (l - p 0 2 )\ 


(36) 


n 2 > S<#>2 R 2 
Po ^ r 2 Ri 

Sy 


(37) 



condition is obtained. Under this condition, the proposed 
estimator is more efficient than simple estimates [8]. 

4.5. Comparison of the Proposed Estimators with Each 
other 

In the case of 

A[R 1 2 Sf+S y 1 (l-pf)<A[R j 2 Sf+S y 2 (l- Po 2 )]; j = 2,3 (38) 

D 2p 2 2 C 2 

I?! V - Rj 2 S^ 2 < 0 

S/(R 1 2 -R J 2 )<0 

l*i I < |*; | (39) 

t 1 estimate is more effective than t 2 and t 3 . 

4.6. Comparison of the Proposed Singh, Chaucer, Sawan, 
Smarandache Estimates and Naik and Gupta 
Estimations 

With the proposed estimation, comparison of Naik and 
Gupta estimate is achieved under the condition of 

MSE (y NG ) - MSE(t t ) >0 ; i = 1, 2, 3 (40) 



From Equation (11) and (18); Under the condition 

A[(Rkci 2 S x 2 + S y 2 ( 1 - P 2 )] < A(S y 2 - 2RpSyS x + R 2 S 2 ) 

{Rk a 2 ~ R 2 )S 2 < p 2 Sy 2 - 2 RpSyS x ) 

{Rkpi 2 - R 2 )S 2 < pSy(pSy - 2RS X ) (32) 

Kadilar-(,'ingi estimation is more effective. 

For example, consider the case i = 1 
If r kc1 = r = l 
Under the condition 

0 < pSy(pS y — 2 RS X ) 2 RS X < pSy 



Po 2 > yj 
•V 



R 2 



■ Rep 2 + 2 RqPo t- 



( 41 ) 



Under this condition, mean square error of the proposed 
estimator is less and therefore more effective [8], 



4. 7. Comparison of Naik and Gupta Estimate and Simple 
Estimation 



ASy 2 > A(R ng 2 Sq 2 - 2 R NG Sy$ + Sy 2 ) (42) 



2 PoSyS^ > RncSq 2 (43) 



Po > 



R NG s cf, 
2 Sy 



(44) 



p>2 C f (33) 

L y 

Kadilar-(, ingi estimation is more effective. 

4.4. Comparison of the Proposed Singh, Chaucer, Sawan, 
Smarandache Estimations and Simple Estimate 

Comparison of the proposed estimators and simple 
estimates and variance of the sample mean was examined 
in simple random sampling as follows; 



condition is obtained. 

Under this condition, we can conclude that mean square 
error of Naik and Gupta estimation is smaller and therefore 
more effective. 

5. Application 

School ranking of 18 districts in Trabzon province in the 
2011-2012 academic years was used (Trabzon MEM, 2011- 
2012) [9]. The data from 18 districts of Trabzon province 
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included the number of students and teachers from 
secondary schools in the 2011-2012 academic years. 
Information regarding the number of teachers and students 
was obtained from Trabzon provincial Directorate of 
National Education. 

Main variable (y) is taken as the number of teachers and 
auxiliary variable (x) is taken as the number of students. 



The relationship between the number of teachers and the 
number of students can be demonstrated by a linear 
equation (y = bx) passing through the point of origin and 
the related variables can be estimated using auxiliary 
variables. Mass unit is the 111 secondary schools in 
Trabzon province. Data regarding the mass is given in 
Table 5.1. 



Table 5.1. Mass units regarding the number of teachers (y) and the number of students (x) variables. 











N = 111 


F = 29.28 


X = 395.06 


p = 0.939 


Y 

R = ■= = 0.07 
X 


S y 2 = 651.312 


S y = 25.521 


Sy 

C v = ^- = 0.872 
y Y 


S x 2 = 159948.937 


S x = 399.936 


s x 

C x = -^r= 1.012 
* X 


r 

ii 

*l|4 /5 

ii 

o 

CO 

NJ 

sD 



Samples are obtained from a mass via simple random 
sampling without replacement using SPSS-20 program. 
First, sample size is calculated. 



d 2 (N-l)+t 2 PQ V 

can be used in the estimation of sample size [10]. 

In this equation, P ratio of units having the same features 
is (Q = 1 — P). If we do not have information about units 
ratio which has the same feature, the largest sample that can 
be selected for P = Q = 0.5 could be reached. lTowever, 
unit ratio having specific features is calculated as the ratio 



of the schools in the city, where d is the amount of 
tolerance. Considering the above equation (45), sample size 
reduces when the the amount of tolerance increases. The 
amount of tolerance (d) depends on the researchers. In the 
present application, if d = 0.09 sample size is obtained 
as 



n = 



1 , 96 2 * 111 * 0 . 81 * 0.19 
0 . 09 2 * 1 10 + 1 . 96 2 * 0 . 8 1 * 0.19 



= 45 



(46) 



The statistics and the data regarding the schools location 
in the city are given in Table 5.2 and Table 5.3 respectively. 



Table 5.2. Sample units of the number of teachers (y) and the number of students (x) variable. 



n = 45 


y = 32.02 


X = 413. 82 


r = 0.955 


1? = — = 0.077 

X 


S y 2 = 875.522 


S y = 29.589 


n 

/ = -= 0.405 


S 2 = 183282.5586 


S x = 428.115 


1-/ 

A = - = 0.013 


b = = 0.066 






n 


S X 



Table 5.3. Mass and sample data required for unit ratio having certain features. 











. . EgB 






mm 











Table 5.4. Mean Ratio Estimators and Mean Square Error 



Mean Estimators 


Mean Values 


Mean Square Error 


Simple Estimation (y) 


32.02 


8.467 


Conventional Proposional Estimation (y 0 ) 


30.568 


2.168 


Singh-Tailor Estimation (ysr) 


30.572 


2.163 


Kadilar-Qmgi- 1 Estimation (y^i) 


29.386 


11.190 


Kadilar-Cingi-3 Estimation (y K Q 3 ) 


29.389 


12.388 


Kadilar-Cingi Estimation ( y prl ) 


29.389 


12.357 


Kadilar-^ingi Estimation (y pr2 ) 


29.389 


12.082 


Kadilar-^ingi Estimation (y pr3 ) 


29.390 


13.880 



According to the mean ratio estimators, the average province was 30. In addition, Singh- Tailor estimation 

number of teachers per secondary schools in Trabzon (y ST ) showed the highest sensitivity. 
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Table 5.5. Ratio Estimators having certain features and Mean Square Error. 



Estimators of units ratio having certain features 


Estimates 


Mean square error 


Naik-Gupta Estimate ( y NG ) 


29.949 


9.058 


Singh-Chaucan-Sawan-Smarandache 1 Estimate(t 1 ) 


29.216 


10.705 


Singh-Chaucan-Sawan-Smarandache 2 Estimate(t 2 ) 


29.940 


9.10 


Singh-Chaucan-Sawan-Smarandache 3 Estimate(t 3 ) 


29.619 


9.713 



According to ratio estimators having certain units feature, 
the average number of teacher per secondary schools in 
Trabzon was 29. Naik-Gupta Estimate (y NG ) estimator was 
the most efficient estimator with the minimum mean square 
error. 

Comparison of Singh and Tailor Estimataion with 
Conventional Proportional Estimation 

MSE(y ST ) = A(S y 2 + R 2 S x 2 w(w - 20)) = 2.163 

where w = =3- = 0.998 and 0 = — = 0.809 
X+p c x 

To ensure inequality MSE(y ST ) < MSE(y 0 ) w — 1 < 0 
and w+1 — 20>0 or w — 1 > 0 and w+1 — 20<O 
condition should be fulfilled. Thus, since w — 1 = 0.998 — 
1 < 0 and w+ 1 - 20 = 0.998 + 1 - 2 * 0.809 = 
1.967 > 0 condition is satisfied, 

MSE(y ST ) = 2.163 < MSE(y 0 ) = 2.168 

mean square error of Singh and Tailor Estimation is smaller 
and therefore more effective. 

Comparison of Singh and Tailor Estimation and Simple 
Estimation 

The variance of simple estimation is V(y) = AS y 2 = 
8.467 

MSE(y ST ) < MSE(y) inequality is valid when the 
condition — = = 1.234 < 2 is satisfied. Since this 

0 0.809 

condition is fulfilled, the Tailor Singh estimation is more 
effective than simple classical estimation. 

Comparison of Kadilar-ffngi-1 Estimation and 
Proportional Estimation 

Under the condition MSE^k^) = A.[(R 2 S X 2 + 

Sy21— p2=11.190 p>iCxCy, Kadilar 

^mgil estimation is more effective than conventional 
proportional estimate. However, condition is not satisfied, 
so conventional proportional estimate is more efficient. 

Comparison of Naik-Gupta Estimation and Simple 
Estimation 

Mean square error is found as MSE(y NG ) = 
X(R ng % 2 — 2 R NG S y( p +5 y 2 ) = 9.058. When compared 

with simple estimates, if the condition p 0 > * is 

satisfied, Naik and Gupta Estimation is smaller and 
therefore more effective than simple estimation. However, 

since p 0 = 0.216 < NG ^ condition is not satisfied, error 

2 Sy 

of simple estimate is smaller and therefore more effective. 

Comparison of the Proposed tq Estimation and Simple 
Estimation 

To ensure V(y) — MSEfj) > 0 inequality, p a 2 > 

s 2 

2 condition must be fulfilled. 

q Z L 

b y 



For this condition i — 1; 



p 0 2 = 0.216 2 = 0.046 < 



0.155 * (36.148) 2 
651.312 



0.311 



s ^ 

Since p a 2 > -^jR^ ondition is not satified, error of 

s y 

simple estimation is smaller and therefore more effective. 
For i = 2; 



p 0 2 = 0.216 2 = 0.046 < %R 2 * = 

Sv 



0.155*(22.593) 2 



651.312 

0121’dir. Therefore, simple estimation is more effective. 
In the same way, simple estimate is more effective for 
i — 3 case. 

Comparison of the proposed tq estimation and t 2 
estimation 



Since /?,=-= 36.15 > R 2 = — = 22.59, 

1 p 2 p+c p 

MSE(t 2 ) = 9.10 < MSE(tf) = 10.71 error of t 2 
estimation is smaller and therefore more sensitive. 



6. Conclusion 

When B = —Rgpi in Kadilar-t^mgi (2004) and B = 
—Rpri in Kadilar-Qingi (2006) estimation, mean square 
error of Kadilar-t^ingi (2004) and Kadilar-^ingi (2006) 
estimates are equal to that of simple estimation. 

When comparing Singh and Tailor Estimator with the 
conventional proportional estimators, under the condition 
w — 1 < 0 and w + 1 — 26 > 0 or w — 1 > 0 and 
w + 1 — 2.6 < 0 Singh-Tailor estimation is more effective 
than conventional proportional estimation. 

In our numerical example, since 
w — 1 = 0.998 — 1 < 0 and w + 1 — 26 — 0.998 + 1 — 
2 * 0.809 = 1.967 > 0 condition is satisfied, MSE(y ST ) < 
MSE(y 0 ) inequality is obtained. 

In the comparison of Singh and Tailor Estimators with 
the simple estimate, j < 2condition is satisfied and under 
this condition, Singh-Tailor Estimator is more effective 

, . , . T ... w 0.998 

than simple estimator. In our applications, — = = 

1.234 < 2 condition is fulfilled and Singh-Tailor Estimator 
was found to be more sensitive. 

In the ratio of units having certain features, estimators 
proposed by Singh et al., tq , t 2 and t 3 estimators were 
compared and R r < Rj (j = 2,3) condition was obtained. 
Under this condition, the proposed Iq estimation is more 
effective than t 2 and t 3 estimations. In our application, 

R 1 36.15 > R 2 — — — — 22.59 inequality is 

P P+Cp 

valid and expressed as MSE(t 2 ) < MSE(tf). 

Considering the numerical examples, according to the 
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estimators proposed in the literature, the average number of 
teachers per secondary schools in Trabzon province in 
2011-2012 academic years ranged from 29 to 32. Singh- 
Tailor estimation ( y ST ), is the mean ratio estimator having 
the smallest root mean square error (2,163) whereas the 
estimator proposed by Kadilar and £tngt ( y pr3 ) has the 
greatest mean square error (13,880). Among the estimators 
having certain features, Naik-Gupta estimator (9.058) has 
the smallest mean square error. 

The numerical results obtained are valid only for this 
application. 
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